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ABSTRACT 

Aims. The statistical properties of planets in binaries were investigated. Any difference to planets orbiting single stars can shed light on 
the formation and evolution of planetary systems. As planets were found around components of binaries with very different separation 
and mass ratio, it is particularly important to study the characteristics of planets as a function of the effective gravitational influence 
of the companion. 

Methods. A compilation of planets in binary systems was made; a search for companions orbiting stars recently shown to host 
planets was performed, resulting in the addition of two further binary planet hosts (HD 20782 and HD 109749). The probable original 
properties of the three binary planet hosts with white dwarfs companions were also investigated. Using this updated sample of planets 
in binaries we performed a statistical analysis of the distributions of planet mass, period, and eccentricity, fraction of multiplanet 
systems, and stellar metallicity for planets orbiting components of tight and wide binaries and single stars. 

Results. The only highly significant difference revealed by our analysis concerns the mass distribution of short-period planets. Massive 
planets in short period orbits are found in most cases around the components of rather tight binaries. The properties of exoplanets 
orbiting the components of wide binaries are compatible with those of planets orbiting single stars, except for a possible greater 
abundance of high-eccentricity planets. The previously suggested lack of massive planets with P > 100 days in binaries is not 
confirmed. 

Conclusions. We conclude that the presence of a stellar companion with separation smaller than 100-300 AU is able to modify the 
formation and/or migration and/or the dynamical evolution history of giant planets while wide companions play a more limited role. 

Key words. (Stars:) planetary systems - (Stars:) binaries: visual - (Stars:) individual: HD 20782, HD 20781, HD 109749, GL86, 
GL86B 



1. Introduction 

The study of the frequency of planets in binary systems' and of 
the properties of these planets and of the binary systems them- 
selves is very important to improve our knowledge of planet for- 
mation and evolution. On one hand, the frequency of planets in 
binary systems has a strong effect on the global frequency of 
planets, more than half of solar type stars being in binary or mul- 
tiple systems (Duquennoy & Mayor 1991). On the other hand, 
the properties of planets in binaries, and any difference to those 
of the planets orbiting single stars, would shed light on the ef- 
fects caused by the presence of the companions. 

The search for planets in binaries can follow two comple- 
mentary approaches. The first one is to perform dedicated sur- 
veys looking for planets in binary systems. Several programs 
are currently in progress focusing on different types of binaries. 
The SARG planet search has been studying about 50 wide pairs 
with similar components (Desidera et al. 2006 ) since late 2000. 
A similar study was recently started by the Geneva group us- 
ing HARPS. About 100 spectroscopic binaries with orbital peri- 
ods between 1.5 and 100 yr were monitored by Eggenberger et 
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' Some of the systems considered in this study are actually hierar- 
chical triple systems. For simplicity, we refer to the whole sample as 
binaries. 



al. r2006a}, while Konacki J2005a> is monitoring about 450 stars 
in binaries. 

The second approach is to study the binarity of the hosts of 
planets discovered in general surveys, which include many bi- 
nary stars in their list in spite of some selection biases against 
them. Dedicated adaptive optics observations are required to re- 
veal faint companions close to the star (e.g. Patience et al.j2005| 
Chauvin et al. 2006 Mugrauer et al. 2006a see Table^for fur- 
ther references). However, the search for brighter companions 
at large separation can be performed by checking astronomical 
catalogs for nearby, comoving objects. 

The first analysis of the properties of planets in binaries re- 
vealed differences with respect to those of planets orbiting sin- 
gle stars. Zucker & Mazeh (2002i showed some difference in 
the period-mass relation. A more complete statistical analysis by 
Eggenberger et al. (2004), hereafter E04, confirms that massive 
planets in close orbits are mostly found in binaries. They also 
propose other possible peculiar characteristics of planets in bi- 
naries: low eccentricities for orbital periods shorter than 40 days 
and a lack of massive planets with periods longer than 100 days 
in binaries. The estimated significance of these latter features is 
less than 97%, therefore requiring larger samples for confirma- 
tion. 

The continuing observational efforts of planet searches and 
dedicated follow-up to study the characteristics of planet hosts 
lead to a significant increase of the total number of planets in 
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binaries. We now know more than 40 planets in binaries or mul- 
tiple systems, an increase by more than a factor of two with re- 
spect to the first compilation of the planets in binary systems 
assembled by E04 (19 planets in 15 systems^). 

The recent publication of the Catalog of Nearby Exoplanets 
by Butler et al. (|^06 ), hereafter B06, provides an updated and 
homogeneous assembly of planet and stellar properties (exclud- 
ing binarity). This work is a useful starting point for a statistical 
comparison of the properties of planets orbiting the components 
of multiple systems and single stars. In this paper we present 
such a discussion. 

The outline of the paper is the following. In Sect. |2 the 
samples of planets in binaries and in single stars are defined; 
in Sect. |5] the statistical analysis of the properties of planets 
in binaries and single stars is presented. In Sect. 0] the results 
and their implications are discussed. The appendices lAICI ad- 
dress specific issues relevant for individual cases: in Appendix 
|X]we present the results of the search for companions of stars 
recently shown to host planets, resulting in two further binaries 
(HD 20782 and HD 109749), whose properties are discussed is 
detail; in Appendix 151 we discuss the probable original proper- 
ties (mass ratio, separation) of the three systems including white 
dwarf (hereafter WD) companions; in Appendix O we provide 
details about individual objects. 



2. An updated compilation of planets in binaries 

We assembled the properties of the planets in multiple systems 
and of the multiple system themselves (Table 0. For the planet 
hosts not included in the latest study of multiplicity of planet 
hosts (Raghavan et al. 2006 hereafter R06), we performed a 
search in existing catalogs of binaries. ^. Details are presented 
in AppendixIXI 

The planetary properties, as well as mass and metallicity of 
the planet hosts, were taken from B06. The binary parameters 
were taken from the references listed in the Table caption. As 
we are interested in evaluating the actual dynamical effects of 
the companions, we listed individual masses, derived using the 
mass-luminosity relations by Delfosse et al. ( 2000) or taken from 
the works cited in the caption. In a few cases, the complete bi- 
nary orbit is available, making feasible a more in depth analysis 
of the possible formation and evolution of the planet (see e.g. 
Thebault et al. 2004 for y Cep). The last column lists the crit- 
ical semiaxis for dynamical stability of planetary orbit derived 
using Eq. 1 of Holman & Weigert ( 1999 1. When the full binary 
orbit is not available, this was calculated deriving the semimajor 
axis from the projected separation and adopting an eccentric- 
ity of 0.35 (Fischer & Marcy 1992 , Duquennoy & Mayor 1991 ). 
Remarks on individual objects are reported in AppendixICl while 
the properties of the systems containing a white dwarf are dis- 
cussed in AppendixlEI 

The comparison sample of planets in single stars is that of 
B06, excluding the objects in Table [2 The mass limit adopted 
by B06 (msin/ < 24 Mj) is different to those adopted in R06. 

^ The radial velocity variations of one of the stars listed in E04 (HD 
219542B) was shown to be due to the activity jitter of the star and not 
to a planet (Desidera et al. 2004^. 

^ Some planets have been announced after the submission of this pa- 
per and they are not considered in this study. One of them is orbiting 
a component of a wide binary (ADS 16402B; Bakos et al. 2006b I. The 
binary separation is 1550 AU. The two transiting planets discovered by 
Collier Cameron et al. 1 2006 1 have companions candidates discovered 
using adaptive-optics imaging, whose physical association to the planet 
hosts needs confirmation. 
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Fig. 1. Mass ratio vs semimajor axis of the binary orbit for stars 
with planets in binary systems. Open circles represent the pairs 
for which binary orbit is available, open squares the pairs for 
which only the binary separation is available. 



One 'super-planet' in the range m sin ; 13 -24 Mj was found in 
a binary, the companion orbiting HD 41004B. 

3. Statistical analysis 

Fig-IDshows the mass ratio vs semimajor axis for stars with plan- 
ets in multiple systems. For hierarchical triple systems in which 
the planet orbits the isolated companion, the masses of the binary 
companions to the planet host are summed. It appears that plan- 
ets might exist in binaries with very different properties. In some 
cases (e.g. very low mass companions at a projected separation 
larger than 1000 AU) the dynamical effects of the companion 
on the formation and evolution of the planetary system might be 
very limited, while in the cases of very tight binaries the pres- 
ence of the planet represents a challenge for the current models 
of planet formation (Hatzes & Wuchterl 2005 1. 

A simple look at Fig.^suggests a possible depopulated re- 
gion for the binary semimajor axis between about 20 to 100 AU. 
A proper evaluation of its reality requires a detailed study of the 
selection effects of the original samples of radial velocity sur- 
veys and the discovery of stellar companions of planet hosts, 
which is postponed for a future study. 

The dynamical effects of the companions on the circumstel- 
lar region of the planet hosts are certainly very different for 
the binaries in our sample. To consider the effects of dynami- 
cal perturbation by the stellar companion(s) we used the critical 
semiaxis for dynamical stability of the planet acii (Holman & 
Wiegert'1999 l. This allows us to take into account the effects 
of both the separation and the mass ratio. The critical semiaxis 
for dynamical stability represents the limit for stable planetary 
orbits around the planet hosts (calculated for coplanar circular 
orbits). It is also similar to the radius of tidal truncation of the 
circumstellar disk (Pichardo et al. 120051 Pfahl & Muterspaugh 
2006 1. However, the area in which relative impact velocities be- 
tween planetesimals is expected to allow planet formation might 
be significantly smaller than flc„-, (Thebault et al. .2006i> . 

The critical semimajor axis flcnr was used to divide the sam- 
ple according to the relevance of the dynamical effects. We de- 
fine as 'tight' binaries those with flf„, < 75 AU and 'wide' bi- 
naries those with Ocrit > 75 AU. The limit corresponds to a pro- 
jected separation of about 200-300 AU depending on the mass 
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Fig. 2. Projected mass vs orbital period of extrasolar planets. 
Open circles: single stars; filled circles: binary stars. The size 
of the symbol is proportional to the critical semimajor axis for 
dynamical stability (larger symbols refer to the tighter binaries). 



ratio. The choice was fixed as a compromise, to obtain a sample 
of binaries wide enough to guarantee limited dynamical pertur- 
bations for the portion of the disk on which the formation of gi- 
ant planets is expected to occur and to allow the inclusion of an 
adequate number of objects in the 'tight' binary sample. When 
relevant for the statistical analysis, the effect of changing this 
limit is considered in the following discussion. 

The statistical comparison to test the hypothesis that the pa- 
rameters of planets (mass, period, eccentricity) in tight and wide 
binaries and in single stars can be drawn from the same par- 
ent distribution was estimated using the Kolmogorov-Smirnov 
(hereafter KS) test and the Mann-Whitney U (hereafter MWU) 
test (Babu & Feigelson 1996 1. We also used the hypergeometric 
distribution (see E04), that yields the probability of obtaining by 
chance a number of cases from a subsample of given size. The 
statistical analysis was performed including only planets with 
radial velocity (hereafter RV) semi-amplitude A' > 15 m/s, to 
exclude planets with uncertain orbital parameters and to ensure 
a more homogeneous detectability. This choice has little effect 
on the results. 



To test this hypothesis statistically, we selected from Tabled 
and B06 the planets with f < 40 days (about in the middle of 
the so-called 'period valley' on the period distribution of plan- 
ets, see Udry et al. 2003 and Fig. |5}- We performed a KS test 
comparing the mass distribution of planets with f < 40 days 
orbiting single stars, binaries with Gcrii < 75 AU and binaries 
with Ocrii > 75 AU (34, 5, and 6 planets respectively) (Fig. |3}- 
The hypothesis that the mass of close-in planets in tight bina- 
ries can be drawn from the same parent distribution of close-in 
planets orbiting members of wide binaries or single stars can be 
rejected with 98.2% and 99.5%. The MWU test yields an even 
higher significance (> 99.9%). The mass distribution of close-in 
planets orbiting single stars and wide binaries is not statistically 
different. 

We checked the effects of the assumptions made on our anal- 
ysis on the resulting significance. The significance of the differ- 
ence of the mass distribution between planets in tight binaries 
and single stars is larger than 99.9% for Ucrit^SO AU, it is about 
97.5% for £3^.^,7— 100 AU and becomes smaller moving the limit 
beyond 100 AU. The exclusion of super-planets more massive 
than 10 My has some eff'ect but the significance remains high 
(> 99%) for both tests. Changing the period limit of close-in 
planets in the range 20-200 days keeps the significance always 
larger than 99%. The RV semi-amplitude selection limit plays 
little role. 

For some of the close-in planets with msin/ > 1.5 My or- 
biting stars classified as single there are indications of the pres- 
ence of significant dynamical perturbations. HD 118203 b, HD 
68988 b and HIP 14810 b have an unusually high eccentricity 
(e - 0.31, 0.14, and 0.15 respectively) for their orbital period 
(~ 6 days). A second planet is orbiting HIP 14810 (B06), while 
in the other two cases linear trends in the residuals from the short 
period orbit indicate the presence of a further companion in a 
long period orbit. Their mass and period will be revealed by the 
continuation of the observations, allowing a proper classification 
of their hosts as binary systems or multi-planet hosts. A candi- 
date companion close to HD 162020 was reported by Chauvin et 
al. ( 2006 1, whose physical association remains to be established. 
For HD 73256 b there is currently no indication of additional 
companions. 



3. 1 . The mass distribution of close-in planets 

Zucker & Mazeh (2002 1 showed that planets in binaries follow a 
different mass period relation than those orbiting single stars. In 
particular, high-mass, short-period planets appear to be present 
only in binary systems. Such a correlation was confirmed by £04 
but further inclusion of new systems makes it weaker (Mugrauer 
et al. 2005 1. The addition of further low-mass hot-Jupiters in bi- 
naries from our work consolidates this latter trend. 

The possibility that the Zucker & Mazeh correlation was 
simply due to the small number of planet in binaries known 
at that time can be considered. However, an alternative inter- 
pretation appears more convincing when one takes into account 
the binary properties of the massive hot-Jupiter hosts: these 
are mostly relatively tight binaries (Fig. |2jl. Indeed, Zucker & 
Mazeh (2002 1 considered in their statistical analysis only bina- 
ries whose projected separation is smaller than 1000 AU, then 
excluding very wide binaries included in increasingly larger 
number in later studies. 



3.2. The mass distribution of planets in wide orbits 

Fig- S shows the distribution of msin; of planets with f > 40 
days. The lack of massive planets with P > 100 days in bina- 
ries noted by E04 is no longer valid: three planets with m sin / > 
5.0 My are now in the sample. All of them are orbiting the com- 
ponents of very wide binaries (projected separation larger than 
800 AU). The probabihty of having 3 planets with P > 100 days 
and msin/ > 5.0 My derived using the hypergeometric distri- 
bution is 17% when considering only wide binaries and 7.4% 
when considering all the binaries together However, this prob- 
ability becomes larger than 10% when considering a different 
mass limit (4-6 My) and KS and MWU tests do not reveal 
statistically significant differences. 

The mass distribution of long-period planets around tight bi- 
naries (7 planets) is not statistically significantly different with 
respect to that of planets in single stars and wide binaries (83 
and 19 planets respectively). 
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Fig. 3. Distribution of the projected mass msin/ for the plan- 
ets with period shorter than 40 days and K < 15 m/s orbit- 
ing single stars (continuous line), planets orbiting binaries with 
ttcrit < 75 AU (dashed line), and planets orbiting binaries with 
dcrit > 75 AU (dotted line). Upper panel: histogram; lower panel: 
cumulative distribution. 



Fig. 4. Distribution of the projected mass m sin / for the planets 
with period longer than 40 days orbiting single stars (continous 
line), planets orbiting binaries with acru < 75 AU (dashed line), 
and planets orbiting binaries with Qcni > 75 AU (dotted line). 
Upper panel: histogram; lower panel: cumulative distribution. 



3.3. The period distribution 

The period distribution of short period planets shows a pile-up 
for periods around 3 days (see e.g. B06). Only one of the five 
close-in planets in tight binaries has a period close to 3 days 
(t Boo b). Two have a very short period (HD41004Bb and HD 
189733b have the shortest and the third shortest period in the 
B06 catalog), while GL86b and HD 195019b have longer peri- 
ods (15 and 18 days respectively). The differences are not statis- 
tically significant using KS and MWU tests. 

When considering the period distribution of planets with P > 
40 days, there are hints that the period distribution of planets 
in tight binaries is different, with a lack of planets with periods 
longer than 1000 days. However, the KS and MWU tests do not 
show statistically significant differences. 

When considering the whole period range (Fig. |5}, the dif- 
ference between the period distribution of planets around single 
stars and tight binaries (117 and 12 planets respectively) is more 
significant (94.3%) for the MWU test while it remains not signif- 
icant according to the KS (50%). Similar values of significance 
are derived for the comparison between the period distribution 
of planets in tight vs planets in wide binaries (25 planets). The 
discrepant level of significance between the two tests can be un- 



derstood considering that the period distribution of planets in 
tight binaries appears to be shifted with respect to those orbiting 
single stars for the whole period range, but without very large 
deviations. The maximum difference of the cumulative distribu- 
tions, the estimator used by the KS test, is then not significant, 
while the global shift is marginally significant by the MWU-test, 
that is sensitive instead to differences of the median of the distri- 
butions. 

As the significance is not large and the different surveys have 
their own period sensitivity and biases concerning binarity (see 
Sect, \3.1i . we conclude that the evidence of the different period 
distribution for planets orbiting single stars and components of 
tight binaries needs confirmation. 

The period distribution of planets orbiting single stars and 
components of wide binaries are instead remarkably similar con- 
sidering both separate period ranges and for the whole period 
range. 

3.4. The eccentricity distribution 

The possible link between high planet eccentricity and binarity 
has been considered since the discovery of the planet orbiting 16 
Cyg B (Holman et al. l 19971 Mazeh et al. il997< . It is thus relevant 
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Fig. 5. Distribution of log(period) for the planets orbiting sin- 
gle stars (continous line), planets orbiting binaries with Oeri, < 
75 AU (dashed line), and planets orbiting binaries with acii > 
75 AU (dotted line). Upper panel: histogram; lower panel cumu- 
lative distribution. 



to consider the eccentricity distribution for planets orbiting tight 
and wide binaries and single stars. 

Fig |6l shows the period-eccentricity diagram. The only two 
planets with e > 0.80 (HD 80606b and HD 20782b) are in binai-y 
systems. However, planets with e up to 0.75 appear to be present 
also around stars currently known as single. 

The eccentricity distribution of planets in tight binaries with 
periods longer than 40 days (Fig.0 is not significantly different 
to those orbiting single stars'*. 

The comparison for planets in wide binaries indicates no sig- 
nificant differences according to the KS test (74%), while the 
MWU test gives a significance of 91.2% that the median of the 
two samples are different. 

E04 noted that all the five planets in binaries in their sample 
with periods shorter than 40 days have an unusually low eccen- 
tricity (e < 0.05). They estimated a chance probability of 3.8% 
using the hypergeometric distribution. While a tendency to a low 
eccentricity of short-period planets in binaries seems to be still 
present (4 planets of 5 and 4 of 6 with e < 0.05, for the tight and 
wide binaries respectively, compared with 17 of 34 for planets 
in single stars) a similar analysis in our sample (taking together 
wide and tight binaries) yields a probability of 12%, therefore 
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The large eccentricity of HD41004Ba is uncertain (e = 0.74±0.20). 



Fig. 6. Eccentricity vs orbital period for planets in binaries (filled 
circles) and orbiting single stars (empty circles). Different sizes 
of filled circles refer to different periastron of the binary orbit 
(larger sizes: closer orbits). 



not confirming the E04 result at a highly significant level. KS 
and MWU tests also do not suggest statistically significant dif- 
ferences. Selecting only planets with periods between 10 and 40 
days (to exclude objects whose orbits are circularizated by tidal 
effects) does not increase the significance of a lower eccentricity 
for short-period planets in binaries. 



3.5. Multi-planet systems 

The occurrence of systems with more than one planet around the 
components of wide binaries (4 of 23, 17.3 + 12.3%) is similar 
with respect to that of planets orbiting single stars (15 of 112, 
13.3 + 4.7%). On the other hand ,no multiple planets have been 
yet discovered around the components of tight binaries (13 ob- 
jects). The closest binary with a known multi-planet system is v 
And (projected separation 750 AU). However, the small number 
of tight binaries in the sample makes the lack of multi -planet sys- 
tems not highly significant (probability of 15% of occurring by 
chance using the hypergeometric distribution). The similar frac- 
tion of multiplanet systems around single stars and wide com- 
panions is a further indication that the presence of a wide com- 
panion does not alter too much the process of planet formation. 



3.6. Metallicity of planets in binary systems 

The study of the metallicity of planet host in binaries can shed 
light on possible different planet formation mechanisms (see 
Sect. 1^. The close-binary planet-hosts appears slightly more 
metal poor than single stars, while the wide-binary planet-hosts 
are slightly more metal rich. The comparison of cumulative dis- 
tributions (Fig.|8} using KS and MWU tests indicates probabili- 
ties of 12.1% and 4.6% respectively that the metallicity of close 
and wide binary planet hosts can be drawn from the same par- 
ent population. The differences with respect to single stars are 
less significant (probability of 23.3% and 17.8% for the close 
binary vs single stars and and wide binary vs single stars respec- 
tively using the MWU test, a smaller significance using KS). 
Considering together all the binaries makes any difference not 
significant. 
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Fig. 7. Cumulative distribution of the planet eccentricity for the 
planets with period longer than 40 days orbiting single stars 
(continous line), planets orbiting binaries with Ocnv smaller than 
75 AU (dashed line), and planets orbiting binaries with Gent 
larger than 75 AU (dotted line). 



Fig. 8. Distribution of the stellar metallicity [Fe/H] for the plan- 
ets orbiting single stars (continous line), planets orbiting binaries 
with Ucrit < 75 AU (dashed line), and planets orbiting binaries 
with fl„„ > 75 AU (dotted line). Upper panel: histogram; Lower 
panel: cumulative distribution. 



While intriguing, the possible differences of the metallicity 
of binary planet hosts should be taken with caution. Beside the 
rather large probability of a chance result, it should be consid- 
ered that the small number of planets in close binaries makes the 
role of individual objects critical (e.g. HD 1 14762, whose com- 
panion might be a brown dwarf or even a low mass star, Cochran 
et al. ll99 1V Furthermore, indirect selection effects might cause 
spurious coiTelations (see Sect. l3.7l i. 

3.7. Selection effects and caveats 

The B06 catalog includes planets discovered by several different 
surveys. Different RV surveys typically have their own biases 
concerning binarity and other properties (e.g stellar mass and 
evolutionary phase, metallicity, etc.). Furthermore, they have 
different precision and time spans, and thus different planet mass 
and period sensitivities. 

All these effects cannot controlled in our study, and we can- 
not exclude that they play a role in producing some of the 
marginal differences we showed (e.g. in metallicity and the pe- 
riod distributions). However, the most significant result, the dif- 
ferent mass distribution of close-in planets in tight binaries, 
should not be affected by these selection effects. In fact, these 



planets are the easiest to discover for any survey because of their 
short periods and large radial velocity amplitude. 

An additional caveat is represented by the uncertainty in the 
upper limit of planetary mass, which is assumed in this study to 
be 24 My following B06, and by the other choices performed in 
the analysis (e.g. period limits, etc.) 

The sample of binaries among planet hosts is certainly not 
complete. A large fraction of the planet hosts were scrutinized 
only for comoving companions at rather large separations (a few 
arcsec). The analysis was performed considering the actual bi- 
nary configuration that in some cases can be significantly differ- 
ent to that at the time of the formation of the planets. In the case 
of binary systems including a white dwarf, the probable original 
configuration can be derived estimating the original mass of the 
(now) white dwarf (see AppendixIBt. but other cases are specu- 
lative (see Appendi? jA.l> . 

4. Discussion 

Exploiting the continuously expanding samples of extrasolar 
planets, we investigated the occurrence of a difference between 
planets in binaries and planets orbiting single stars. The size of 
the sample allowed us to divide the binaries in two sub-groups. 
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'tight' and 'wide' binaries. The separation was fixed at a criti- 
cal semiaxis for dynamical stability equal to 75 AU. This cor- 
responds to separations of about 200-300 AU depending on the 
mass ratio between the components. 

4.1. The properties of planets in tigiit binaries 

The most significant feature we were able to show is that the 
mass distribution of close-in planets around close and wide bina- 
ries is different. Only the close-in planets in tight binaries have 
a mass distribution significantly different to those orbiting sin- 
gle stars. This result implies that the formation and/or migration 
and/or dynamical evolution processes acting in the presence of a 
sufficiently close external perturber are modified with respect to 
single stars. 

The fact that the hosts of massive planets in close orbit are 
typically binaries somewhat resembles the results by Tokovinin 
et al. J2006l i.e. short-period spectroscopic binaries have in most 
cases a further companion. That a third member could act as 
an angular momentum sink resulting in a closer inner binary 
was first suggested by Mazeh & Shaham ( 1979 1. Tokovinin et 
al. |2006 ) favor generation of high orbital eccentricity by Kozai 
resonances, and subsequent tidal circularization, as the mecha- 
nism to explain the close spectroscopic binaries and their link 
with the presence of additional companions. 

While it cannot be excluded that some of the massive close- 
in planets in binaries were moved into their current location by 
Kozai migration, some other massive planets (GL 86b and HD 
195019b) have periastrons of 0.108 and 0.137 AU, probably ex- 
cluding the occurrence of tidal circularization. This calls for a 
more general process to explain the presence of massive close- 
in planets. 

A possible effect on migration and mass accumulation 
timescales in a disk in the presence of a companion was pre- 
dicted by the models by Kley (2000 1. He considers the evolution 
of a Jupiter mass planet embedded in a circumstellar disk in bi- 
nary systems with a separation of 50-100 AU, with eccentricity 
ei,i„ = 0.5, and stellar masses of 1.0 and 0.5 Mq. The presence of 
the companion significantly enhances the growth rate and make 
the migration timescale of the planet shorter As discussed by 
E04, this can qualitatively explain the presence of massive close- 
in planets in binary systems. A step forward here is that we have 
shown that most of the binary systems that host the 'anoma- 
lous' massive, close-in planets are indeed not too different to 
those modeled by Kley J2000I I. while E04 were forced to con- 
sider also binary systems with very wide separation, for which 
the described model is likely not to work. The models by Kley 
( 2000 1 have as an initial condition an already formed planet and 
thus do not include the formation phase of giant planets. As dis- 
cussed by Thebault et al. (|2)06(, a possible way to achieve the 
formation of planets in a swarm of planetesimals perturbed by a 
stellar companion might be type II runaway accretion identified 
by Kortenkamp et al. ( 2001 ). More detailed modeling is required 
to achieve a full understanding of the processes that cause the 
differences of planet properties in tight binaries. 

Additional possible differences of the properties of planets in 
tight binaries (to be confirmed as the statistical significance is not 
very high and selection effects might be at work) are the different 
period distributions (there are no planets with P > 1000, the 
close-in planets probably do not pile-up at period close to three 
days) and the lack of multi-planet systems. 

In a few cases, the lack of long period planets can be under- 
stood as the result of dynamical stability constrains: for y Cep 
and HD 41004A, the critical semimajor axis for dynamical sta- 



bility is just 1.9 and ~3 times the planet semimajor axis. For 
the other cases, a more general mechanism should be at work 
(enhanced migration; Kley 2000 or lack of suitable conditions 
for planet formation in the external parts of the disks, under the 
gravitational influence of the companion). 

The fact that properties of planets in tight binaries are differ- 
ent to those orbiting single stars can be exploited to distinguish 
between different scenarios for their origin. As discussed by e.g. 
Hatzes & Wurchterl J2005> and Thebault et al. (2004i, the pres- 
ence of planets in binaries with a separation of 20 AU or less 
represents a challenge for our current view of planet formation. 
Such problems can be overcome if the binary was initially wider, 
and its orbit was modified through dynamical interactions in star 
clusters after planet formation (Pfahl & Muterspaugh 2006 1. The 
different mass distribution of planets in close binaries implies 
that, if this is the case, the dynamical interactions should not 
only modify the binary orbit but also force a massive planet in 
external orbit to migrate inward, as massive hot Jupiters are rare 
around single stars. This point should be addressed by future 
dynamical simulations. Kozai oscillations and tidal circulariza- 
tion after the dynamical interaction cannot be excluded in some 
cases (large relative inclinations between the planet and binary 
orbits might be the result of the proposed dynamical encoun- 
ters). However, less massive planets originally in external or- 
bits should follow a similar fate and thus we expect that low 
mass planets should also be present in close orbits in tight bi- 
naries, at odds with current observational data. Furthermore, as 
discussed above, Kozai migration cannot explain planets with 
periods longer than a few days and low eccentricities. 

Further clues to the origin of planets in tight binaries can 
be derived by a determination of their frequency. Pfahl & 
Muterspaugh ( 2006 ) predict that the fraction of close binaries^ 
that dynamically acquire giant planets is about 0.1%, with an 
uncertainty of about one order of magnitude. They note that the 
rough number of these planets seems too high to be compatible 
with these predictions. However, the number of tight binaries in- 
cluded in the radial velocity samples is not well known. A more 
detailed analysis is postponed to a forthcoming study. 

The mentioned difficulties of the binary-interaction scenario 
to account for the occurrence of planets in close binaries leave 
open the alternative hypothesis, i.e. giant planets formed in bi- 
naries with small separation at the time of planet formation, pos- 
sibly in a different way than planets around single stars. 

The metallicity distribution of tight binaries with planets can 
give clues to the formation mechanism for these planets. The 
large effect of the stellar metallicity on the frequency of plan- 
etary companions is well established (Fischer & Valenti 120051 
Santos et al. 2004 1 and can be understood in the framework of 
the core-accretion model for the formation of giant planets (Ida 
& Lin 12004. Robinson et al. i200 6i. If these planets formed in 
a difi'erent way than those orbiting single stars and the compo- 
nents of wide binaries, e.g. by disk instability triggered by dy- 
namical perturbations (Boss 2006 1, we should expect a lower 
metallicity for planet hosts in close binaries, as the disk instabil- 
ity mechanism is expected to be rather insensitive to metallicity 
(Boss EUn^ . There is a marginal indication of a larger fraction 
of metal poor stars among tight binaries with planets. However, 
the distribution of metallicity remains on average larger than that 
if the solar vicinity, with a mean and median larger than the so- 



' The definition of close binaries used by Pfahl & Muterspaugh 
(2005J (semimajor axis less than 50 AU) is different to that adopted 
here. 



8 



S. Desidera & M. Barbieri: Properties of planets in binary systems 



lar value. This indicates that high metallicity is a factor favoring 
planet formation also in tight binary systems. 

4.2. The properties of planets in wide binaries 

The mass and period distributions of planets in wide binaries are 
not statistically significant different to those of planets orbiting 
single stars. The lack of long-period massive planets in binaries 
shown by E04 is not confirmed at a highly significant level in our 
study. The fraction of multi-planet systems is similar for planets 
orbiting single stars and components of multiple systems. These 
results indicate that a distant companion (separation >300-500 
AU) does not significantly affect the process of planet forma- 
tion. The only marginally significant difference between planets 
orbiting single stars and components of wide binaries concerns 
the planet eccentricity. 

4.3. Is binarity the cause of the high eccentricities ? 

The eccentricity distribution of planets in wide binaries shows a 
marginal excess of high-eccentricity (e > 0.6) planets. However, 
current results indicate that the high planet eccentricity is not 
confined to planets in binaries, and that the possible differences 
in eccentricity are limited to the range e > 0.5 - 0.6. This indi- 
cates that there are mechanism(s) generating planet eccentricity 
up to 0.4-0.5 that are independent of the binarity of the planet 
host, and are characteristic of formation and evolution of a plane- 
tary system (e.g. disk-planet interactions: Tremaine & Zakamska 
120041 planet-planet scattering: Marzari & Weidenschilling i2002i 
Ford et al. ^03 1. These probably act during or shortly after 
planet formation. Further eccentricity enhancements, possibly 
linked to the presence of a companion, might take place at 
later epochs. In fact, Takeda et al. (I2006> noted that most high- 
eccentricity planets orbit old stars (ages >5 Gyr). Mechanisms 
that require long time scales to modify planetary orbits then 
seem favored. 

A popular mechanism that might produce high eccentricity 
planets in binaries is represented by Kozai oscillation. Kozai os- 
cillation is expected to be at work also for low mass companions 
at large separation (the case of most of the companions of planet 
hosts discovered up to now), provided that a sufficiently long 
time is available and that the planet has a sufficiently long pe- 
riod to prevent the suppression of Kozai oscillation by general 
relativistic eff'ects (Takeda & Rasio 2005 , Wu & Murray 2003 1. 
As an example, a brown dwarf companion such as HD 365 IB 
(mass about O.O6M0, projected separation about 500 AU) might 
be able to induce the eccentricity oscillation for a planet with an 
initial semimajor axis larger than about 3 AU. 

Kozai oscillations are expected to produce large planet ec- 
centricities. Wu & Murray (2003 1 considered in detail the pos- 
sibiUty that the high eccentricity of HD 80606b (e - 0.93) is 
due to Kozai oscillation. This was found to be possible only 
for a relative inclination of the planetary and binary orbit larger 
than 85 deg. In the case of the only other planet with e > 0.8 
(HD20782b), the period of the Kozai eccentricity oscillation de- 
rived from Eq. 2 of Takeda & Rasio ( 2005 1 is much longer than 
the Hubble time. Therefore, the extremely high eccentricities of 
HD20782b is unlikely to be due to the Kozai mechanism (see 
Appendix lA. 1> . unless further companions exist in the system at 
smaller separation. 

An alternative to Kozai oscillations is represented by chaotic 
evolution of planetary orbits induced by dynamical perturba- 
tions (Benest & Gonczi .2003. Marzari et al. i2005<) . This kind of 



evolution might also arise for the dynamical interactions within 
the planetary system, without requiring a stellar companion in a 
wide orbit. 

To further investigate the origin of the high eccentricities, we 
consider the eccentricity distributions of planets residing in mul- 
tiplanet systems and isolated planets (Fig.|5}. There is a lack of 
high eccentricity planets in multi-planet systems (as noted also 
by Takeda et al. 120061. When considering planets with P > 40 
days, the comparison of the two distributions reveals a difference 
that is significant to about a 89.1% confidence level (according 
to the MWU test), similar to the results of the planets in wide 
binaries vs those in single stars. This can be understood consid- 
ering that planetary systems with highly eccentric planets are in 
general more subject to close encounters, and then less stable, 
than a system populated with planets in low eccentricity orbits. 
If the chaotic evolution of a planetary system has as a final re- 
sult a planet on a highly eccentricity orbit, it is probable that that 
planet remains alone. On the other hand, Kozai oscillations are 
more likely to take place in single-planet systems, as they are 
suppressed in the presence of significant mutual planet-planet 
interactions. Therefore, eccentricity enhancements due to Kozai 
oscillation are expected to occur for single-planets systems with 
a companion in a suitable binary configuration. A proper eval- 
uation of the link between high eccentricities, binarity of planet 
hosts, and planet multiplicity requires larger samples to allow 
more statistically significant inferences. 

5. Conclusion 

The only highly significant difference shown in our analysis con- 
cerns the mass distribution of short-period planets in tight bi- 
naries. Massive planets in short period orbit are found in most 
cases around the components of rather tight binaries. Other pos- 
sible peculiar features of planets in tight binaries with respect to 
planets orbiting single stars, such as a lack of long-period plan- 
ets and multiple planets, are of marginal significance and need 
confirmation. 

The properties of exoplanets orbiting the components of 
wide binaries are compatible with those of planets orbiting sin- 
gle stars, except for a possible greater abundance of planets on 
highly eccentric orbits. The previously suggested lack of mas- 
sive planets with P > 100 days in binaries is not confirmed. 

This result indicates that the binary separation plays a role 
in affecting the properties of planets and that the simple clas- 
sification of planets in binaries vs planets orbiting single stars 
is not adequate to describe the dynamical effects. This opens the 
perspective for more detailed modeling of the role of stellar com- 
panions on the formation and evolution of planetary systems. 

Appendix A: Binaries among new planet hosts 

As planet announcements arise quite frequently, existing com- 
pilations quickly become incomplete. Twenty-one stars are not 
considered in the latest study of multiplicity of planet hosts 
(R06). Most of them are new planet hosts while some are stars 
with companions with mass between 13 to 24 My, not consid- 
ered by R06 but included in B06 catalog. 

We searched for companions of these stars in existing as- 
tronomical catalogs such as Hipparcos (ESA I1997> . CCDM 
(Dommaget & Nys 2002 1, WDS (Mason et al. lMJB . and ADS 
(Aitken 1932 1. For the following 16 stars there are no evidence 
of binarity in the literature: HD 4308, HIP 14810, HD 33283, 
HD 66428, HD 81040, HD 86081, HD 99109, HD 102195, HD 
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Fig. 9. Distribution of the planet eccentricity for the single plan- 
ets (continous line), and planets in planetary systems (dotted 
line). Upper panel: histogram; lower panel: cumulative distribu- 
tion. 



107148, HD 109143, HD 118203, GJ 581, HD 137510, HD 
187085, HD 212301, HD 224693. 

HD 33564 and HD 164922 are included in CCDM and WDS 
(the latter also in ADS) as triple systems, but the companions 
display different proper motion to the planet hosts and thus they 
are not physical companions. Two additional objects (HD 20782 
and HD 109749) are binaries. The physical association of the 
companions and their physical properties are discussed in the 
next subsections. Finally, the binarity of HD 189733 is discussed 
in a dedicated paper by Bakos et al. ( I2006a> . 

A.1. HD 20782 

HD 20782 is Usted as a binary in the CCDM catalog (CCDM 
03201-2850). The companion is HD 20781. The very large sep- 
aration, 252 arcsec, that corresponds to 9080 AU at the distance 
of HD 20782 (36 pc) makes a detailed check of physical associ- 
ation mandatory. Table lA. ll lists the stellar parameters of the two 
stars. Hipparcos parallaxes and proper motions, and Nordstrom 
et al. (I2004> radial velocities, are fully compatible within errors 
for the two components. This strongly suggests physical asso- 
ciation. Therefore, HD 20782 and HD 20781 form a very wide 
common proper motion pair. Considering the nominal radial ve- 
locity and proper motion differences, the pair could be bound. 



Parameter 


HD 20782 


HD 20781 


Ket. 


fie, (mas/yr) 


348.88+0.50 349.07±0.78 


1 


fjs (mas/yr) 


-64.82±0.73 - 


67.80 ±1.00 


1 


RV (km/s) 


39.5±0.2 


39.6 ± 0.2 


2 


n (mas) 


27.76+0.88 


27.86±1.23 


1 


R,„i„ (kpc) 


4.95 




2 


Rmax (kpc) 


8.16 




2 


ecc 


0.24 




2 


Zmax (kpc) 


0.08 




2 


V 


7.366 


8.457 


2 


ST 


G3V 


KOV 


1 


Mass(Mo) 


0.90 


0.84 


2 




1.00 




3 




969+0.024 




4 


Age;,o,. (Gyr) 


13.0 




2 




7 r' ' 

'•^-4.1 




3 




9.68 ± 1.76 




4 


AgefjK (Gyr) 


3-6 


6 


5,6 



References: 1 Hipparcos (ESA 1997 i; 2 Nordstrom et al. (20(m; 3 
Valenti & Fi scher Ii2005i) : 4 Takeda et al. (2006j; 5 Jones et al. «20051 : 
6 Gray et al. <200gl 



The observed separation is very large but not extreme for bina- 
ries in the solar neighborhood (Poveda & Allen 2004). Clearly, 
the orbit of the wide binary remains unconstrained. 

We tested the dynamical stability of the system via numer- 
ical integration and we have performed 10 simulation of this 
system following the evolution for 200 Myr (577 binary orbits) 
using a RADAU integrator (Everhart et al. 1985 1. The binary 
orbit was assumed with typical values of ai,j„ = 6000 AU and 
^bin - 0.60, the orbital parameter of the planet was chosen equal 
to the observed values, and inclination was chosen slightly in- 
clined (/ = 5°) with respect to the binary orbital plane. The planet 
results stable over all the simulation, thus we argue that the large 
binary separation ensures dynamical stability of the planet in 
spite of its very high eccentricity. 

Only two extrasolar planets have orbits with eccentricities 
larger than 0.8: HD 20782b and HD 80606b. Both their host 
stars are member of wide common proper motion pairs. 

While the link between the high eccentricity and binarity still 
needs statistical confirmation, nevertheless it is interesting to in- 
vestigate the possible ways in which a distant companion might 
have induced extreme planet eccentricities. 

As discussed in Sect. 0] the extremely high eccentricity of 
HD20782b is unlikely to be due to the Kozai mechanism, as 
the period of the eccentricity oscillation is much larger than the 
Hubble time. 

An interesting possibility to explain both the very large sep- 
aration of the binary and the very high eccentricity of the planet 
orbiting the primary is a dynamical encounter of the binary (or 
of an originally higher multiplicity system) within a star cluster 
or with a passing star, that might have perturbed both the binary 
and the planet orbit. According to the simulations by Weinberg et 
al. ( 1987 1, the probability of survival for a binary with an initial 
semimajor axis of about 10000 AU after 4 Gyr is less than 50%. 
Furthermore, the actual semimajor axis is probably larger than 
the observed projected separation by about 30% (Duquennoy 
& Mavor [T99Tl Fischer & Marcv lT992t . further decreasing the 
probability of survival, if this is the case of HD20782. The galac- 
tic orbit of HD 20782 (Table lA.U makes the star on average 
closer to the Galactic center than the Sun, further increasing the 
chance of stellar encounters, because of the higher stellar den- 
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sity. Therefore, at the age of the system (7J_^j Gyr; Valenti & 
Fischer l2005l 9.68 + 1.76; Takeda et al.'2006^, the probabihty 
that the currently observed orbit is not the original one but it was 
significantly modified by dynamical encounters in the Galactic 
disk is rather high. 

The process of destruction of wide binaries is typically the 
result of a large number of weak encounters (see e.g. Portegies 
Zwart et al. .l997j . The changes in the binary orbit caused by one 
or more stellar encounters might alter the orbit of the planet(s) 
orbiting one of the components. A multi-planet system can be 
destabilized, leading to chaotic evolution of the planetary orbits. 
Dedicated dynamical modelling should be performed to check if 
the current configuration of the system might represent a possi- 
ble outcome of a stellar encounter. This is postponed to future 
work. The hypothesis of stellar encounters appears less appeal- 
ing to explain the extremely high eccentricity of HD 80606b, 
as in this case the projected separation of the binary is about 
1200 AU, 7-8 times smaller than for HD 20782. However, the 
real orbit remains unknown. 

The possibility that the high planet eccentricity is unrelated 
to the presence of the stellar companions and is caused by in- 
ternal evolution of the planetary system cannot be excluded. 
Conclusive inferences on the origin of the very high eccentricity 
of HD 20782b and HD 86606b are not possible from the avail- 
able data. Further clues will come from follow-up observations, 
such as the search for further planetary and stellar companions, 
and from the discovery of further planets with very highly ec- 
centric orbit and the study of the binarity of their hosts. 

A.2. HD 109749 

HD 109749 is Usted in Hipparcos and CCDM catalogs as a bi- 
nary (CCDM J12373-4049). The Hipparcos magnitude of the 
secondary transformed to the standard system following Bessell 
(2000 1 yields V - 10.76, about 2.5 mag fainter than the primary 
(V=8.09, Fischer et al. 2006; see Table IX^. Table 1X31 shows 
the relative astrometry for the two components. The constancy 
of the projected separation and position angle, coupled with the 
rather large proper motion of HD 1 09749, makes the physical as- 
sociation between the two objects very probable. The projected 
separation of 8.35 arcsec corresponds to about 500 AU. 

To further check the physical association between the two 
stars, we identified them in 2MASS catalog (Cutri et al. 2003). 
The photometry is shown in Table IA.2I When placed in a 
color-magnitude diagram, both components lies along the same 
isochrone, indicating a common distance and confirming the 
physical association (Fig. IA.1> . The masses derived from the 
isochrone shown in Fig. lA. H are 1.11 and 0.78 Mq. 



Appendix B: White dwarf companions of planet 
hosts 

Three planets are orbiting stars with white dwarfs companions. 
For binaries wide enough to escape common envelope evolu- 
tion, the mass loss from the originally more massive star during 
the RGB and AGE phases determines a widening of the orbit. 
The mass loss phase lasts several Myr, longer than the typical 
binary orbital periods. In these conditions, the widening of the 
binary orbit can be approximated by ae„d = a start M start /M^nd, 
and no eccentricity changes are expected (Lagrange et al. 2006 
and references therein). In this section, we consider the three 
planet-host binaries with a WD, placing constraints on the origi- 



Table A.2. Stellar properties of HD 109749 and HD 109749B 



Parameter HD 109749 HD 109749B Ref. 



^la (mas/yr) 


-157.89±1.41 


-158 


1,2 


fis (mas/yr) 


-5.48+1.25 


-6 


1,2 


n (mas) 


16.94+1.91 




1 


V 


8.09 


10.76 


4,5 


JlMASS 


7.057±0.021 


8.788±0.024 


6 


Him ASS 


6.797±0.031 


8.289±0.051 


6 


K2MASS 


6.678±0.024 


8.123±0.024 


6 


Mv 


4.23+0.25 


6.90+0.25 


1,5 


Mass(Mo) 


1.11 


0.78 


5 


Mass(Mo) 


1.2 




4 



References: 1 Hipparcos (ESA 1997 1; 2 USNO-Bl.O Catalog (Monet 
et al. 2003 1 3 Nordstrom et al. (.2004» : 4 Fischer et al. (l2006l : 5 This 
paper; 6 2MASS (Cutri et al. l2003l 



Table A.3. Relative astrometry of HD 109749 



Epoch 


P 


9 


Ref 




1881 


8.2 


179 


CCDM 




1987.35 


8.35+0.04 


179.76+0.21 


Sinachopoulos 


fT988l 


1991.25 


8.353+0.014 


180.0 


Hipparcos 




1999.27 


8.28 


179.75 


2MASS 
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Fig.A.l. Absolute magnitudes and colors for the components 
of HD 109749 (filled circles). Left panel: vs V-K; right 
panehMy vs J-K. The 4 Gyr, Z=0.030 isochrone by Girardi et 
al. (2002 1 is overplotted. 



nal configuration of the system in terms of binary separation and 
mass ratio. 

B.I. HD 13445 = GL 86 

The companion of G186 was discovered by Els et al. (1200 1> 
and it was classified as a brown dwarf on the basis of the near- 
infrared colors. Mugrauer & Neuhauser (2005 1 showed that the 
companion is instead a white dwarf. Lagrange et al. ( 2006;) fur- 
ther constrained the properties of the system, using archive ra- 
dial velocity and new photometry and position measurements*'. 
However, they also noted some inconsistency between the white 
dwarf cooling age (1.2-1.8 Gyr depending on white dwarf mass 
and atmosphere composition), the age of the star as derived from 
chromospheric activity (about 2-3 Gyr), and the presence of the 
planet. In fact, such a young age implies an original mass of the 
companion that is too large to be compatible with a dynamically 
stable orbit for the planet once the widening of the binary orbit 

* The available observational constraints did not allow Lagrange et 
al. (2006( to derive a unique solution for the binary orbit. The orbit 
listed in TableQis a plausible, representative solution. 
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due to mass loss from the originally more massive star is taken 
into account. 

A picture that might explain the properties of the system is 
the following. Kinematic parameters (U = -102 km/s, V - 
-75 km/s, W - -29 km/s, galactic orbit 7?„„„ = 4.0 kpc, 
R„ax = 9.2 kpc, e = 0.39, Z„„„ = 0.4 kpc; Nordstrom et al. 2004 1 
indicate a very old age^. The kinematic criterium of Bensby et 
al. (2005 1 suggests a probable membership to the thick disk. 
Chemical abundances can also be used to constrain the galac- 
tic population. At the metallicity of GL86 ([Fe/H]~ -0.20), dif- 
ferences between thin and thick stars are fairly small (the typi- 
cal [a/Fe] of thin disk and thick stars is about -h0.05 and -1-0.15 
respectively, Bensby et al. 2 005) . High-resolution abundances 
analysis of several elements were published by Allende Prieto 
et al. (l2004t . Valenti & Fischer (l2005l . Santos et al. (TWAS, 
and Gilli et al. (|7006l. The uncertainty in the chemical abun- 
dances of GL 86 and the discrepancies between different studies 
do not allow a firm classification, but they appear compatible 
with thick disk membership or with thick disk-thin disk transi- 
tion objects^. Thick disk stars with this metallicity have typical 
ages of 8 Gyr according to Bensby et al. (2005 1. This is com- 
patible with the lower limit (68% confidence level) of 8.5 Gyr 
estimated by Takeda et al. ( .200 6 1 from isochrone fitting. If GL 
86 is instead a thin disk star, its kinematic parameters strongly 
supports an old age. 

Matching the cooling age of the white dwarf given by 
Lagrange et al. (2006) with the assumed age of about 8 Gyr for 
the system requires an original mass of about 1 .2 Mq. Older ages 
yield lower initial masses. 

A 1 .2 Mq star ends its life as a ~ 0.54 Mq white dwarf. The 
amount of mass loss implies a widening of the binary orbit of 
about 40-50% of its original separation. To match the semimajor 
axis of 18.4 AU derived by Lagrange et al. i2006> . the semimajor 
axis of the original orbit had to be about 13 AU. 

A fraction of the mass loss by GL86 B during the RGB and 
AGB evolution should have been captured by the planet host. 
A rough estimate based on the BSE binary evolution code by 
Hurley et al. (2002), using the default wind and accretion pa- 
rameters, is 0.017 Mq. The ~ 0.017 Mq accreted mass corre- 
sponds to roughly half of the mass of the convective envelope 
of GL86 (M,, = 0.039 + 0.003; Takeda et al. 2006l. The lack 
of peculiar abundance of neutron capture elements ([Ba/Fe]=- 
0.12; Allende Prieto et al. I2004i) suggests an original mass for 
the AGB star smaller than 1.3 - 1.5 Mq (Busso et al. 1999J , in 
agreement with our estimate of ~ 1.2 Mq. 

Along with mass, angular momentum should also have been 
transferred to the planet host. This was invoked to explain the 
high rotation rates of companions of hot white dwarfs such as 
the barium dwarf 2RE J0357+283 (Jeffries & Smalley .1996). 
The accretion of angular momentum by GL86 would explain its 
activity level typical of a much younger star 

This scenario would be able to explain the available observa- 
tional constraints, including the stability of the planetary orbit. 

' The difference between the kinematic and activity age of GL 86 
was previously noted by Rocha Pinto et al. |20()2|. 

^ The relative abundances of Allende Prieto et al. <2004t 
([O/Fe]=+0.21, [Mg/Fe]=+0.13, [Si/Fe]=+0.08, [Ca/Fe]=-^0.18, 
[Ti/Fe]=-l-0. 11, [Eu/Fe]=+0.30) fit fairly well the expected abundance 
pattern for a thick star of the metallicity of G186, those by Gilli et 
al. 1 2006 1 and Ecuvillon et al. 1 2006 1 are more ambiguous ([0/Fe]=- 
0.01, [Mg/Fe]=+0.21, [Si/Fe]=-0.02, [Ca/Fe]=-0.10, [Ti/Fe]=+0.19). 
Valenti & Fischer (.2005 J find an abundance ratio similar to that 
of Allende Prieto et al. <2004> for two elements ([Si/Fe]=+0.11, 
[Ti/Fe]=+0.10). 



However, the formation of a planet with a companion at about 
13 AU (and with an eccentricity of 0.4) represents a challenge 
for current models (Hatzes & Wuchterl 2005 1. The signature of 
accretion of angular momentum by GL86 is an indication that 
the binary system was in its current status at the time of the mass 
loss by GL86B, about 1.5 Gyr ago. Earlier modifications of the 
binary orbit due to e.g. dynamical interactions within the native 
star cluster (Pfahl & Muterspaugh 2006 1 cannot be excluded. 
The possible effects of the mass and angular momentum accre- 
tion on the mass and orbit of the planet remain to be investigated. 

B.2. HD 27442 = e Ret 

HD 27442= e Ret has a companion at a projected separation of 
about 13 arcsec = 240 AU. According to Chauvin et al. (2006), 
visual and near-infrared photometry are not compatible with 
any main sequence star They then argue that the companion 
is a white dwarf. The mass of the planet host yields a lower 
limit to the original mass of the companion. Valenti & Fischer 
( 2005 ) determine a mass of 1 .49 Mq (the value also adopted by 
B06), with an age of 3.5 Gyr, while Takeda et al. (2006) Hst 
M = 1.48+jj-22 Mq with an age of 2.84+°-^JJ Gyr. A lower limit to 
the WD mass estimated using the BSE code (taking into account 
the high metalHcity of e Ret; [Fe/H]=+0.42) is M ~ 0.60 Mq. 
Theoretical models should be taken with caution at these ex- 
tremely high metallicities, as they fail to explain the WD lumi- 
nosity function of the old super-metal-rich open cluster NGC 
6791 (Bedin et al. 2005 1. The relatively wide separation from 
the primary (13 arcsec) and the moderately bright magnitude 
{V ~ 12.5, WDS) open the perspective for a more detailed char- 
acterization of e Ret B. We can guess that the system had lost 
more than one third of its original mass, implying a substantial 
widening of the binary orbit. The original separation might have 
been about 150 AU, making the influence of the companion not 
negligible. 

B.3. HD 147513 

HD 147413 was classified as a member of the Ursa Major 
moving group by Soderblom & Mayor ( 1983. 1 while King et 
al. (2003) classified the membership as unprobable, on the ba- 
sis of the discrepant kinematic parameters. The activity level 
and the position on the color-magnitude diagram are instead 
fully compatible with membership. Its companion (WD 1620- 
39=H IP 803 00) is a widely studied WD, for which Silvestri 
et al. J200H derived a mass of 0.65 + 0.01 Mq and a cooling 
age of 20 Myr The planet host was claimed to be a barium 
star ([Ba/Fe]=+0.37) by Porto De Mello & Da Silva ( 1997 1 but 
Castro et al. ( 1999 1 showed that the high barium abundance is 
shared by other stars in the UMa moving group (and possibly by 
other young stars). This points toward a primordial abundance 
and not to accretion by the WD companion, unexpected on the 
basis of the large separation (more than 4000 AU). 

We estimated the possible original mass of the (now) WD 
companion using the BSE code (Hurley et al. '2002V Adopting 
the most recent age estimate of the UMa moving group (500 
Myr; King et al. 2003 1 and the cooling age by Silvestri et 
al. (1200 1> leads to an original mass of 3.0 Mq and a WD mass 
of 0.75 Mq, some 0.1 Mq larger than that derived by Silvestri 
et al. (2001V The WD mass and cooling age by Silvestri et 
al. (2001 1 can be matched simultaneously by adopting an ini- 
tial mass of 2.3 Mq and a system age of 1 Gyr. This age is still 
marginally compatible with the photometry and activity level, if 
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the star is not a member of UMa moving group. We conclude 
that the original separation of the binary was about a half of 
the present one, as the system has lost about half of its original 
mass. Considering the current very wide separation, this does 
not change the classification of this planet host as a very wide 
binary. 

Appendix C: Individual objects 

We further discuss individual objects. 

- HD 3651: a brown dwarf companion was recently identi- 
fied by Mugrauer et al. L2006bj . Its mass depends largely 
on the age of the system (20 - 60 My for ages of 1-10 Gyr 
respectively). Takeda et al. (2006) favor a very old age (older 
than 11 Gyr at a 68% confidence level). We then assume 
M = 0.06 Mq. Such a mass coupled with the wide separation 
(480 AU) indicates that the classification of this object within 
the binary companions and not within the 'super-planets' is 
appropriate. The critical semiaxis for this system is uncertain 
as Eq. 1 of Holman & Wiegert (.1999) is not defined for such 
an extreme mass ratio. 

- HD 19994: the visual orbit by Hale J 19941 is very prelimi- 
nary. 

- HD 38529: astrometric acceleration detected by Hipparcos. 
According to RefFert & Quirrenbach (2006.1, this accelera- 
tion is due to the outer 'planet' that has m = 37^]^ My. The 
classification of this star as a two-planet host with a wide 
stellar companion or a single-planet host with a brown dwarf 
and a further wide companion is then ambiguous. 

- HD 41004: the stellar metallicity of HD 41004B is not Usted 
by B06. We assume that of the primary. 

- HD 1 14762: a companion with mass 0.075 Mq was found by 
Patience et al. (2002i. Chauvin et al. (2006 1 report a candi- 
date companion, whose physical association has to be estab- 
lished. 

- T Boo: the very eccentric orbit by Hale ( 119941 e = 0.91) 
implies a periastron of about 20 AU. However, the latest 
position measurement by Patience et al. (2002i indicates a 
larger separation and similar position angle than predicted 
by the preliminary orbit. A redetermination of the binary or- 
bit, taking into account also the long term radial velocity 
trend reported by B06, would be useful. A further L dwarf 
companion candi date at large separation has been reported 
by Pinfield et al. (I2006> . The physical association has to be 
confirmed. 

- HD 178911: hierarchical triple system; a tight pair (masses 
1.10 and 0.79 Mq, semimajor axis 3 AU; Tokovinin et 
al. 2000 1 is at a projected separation of about 640 AU from 
the planet host. 

- HD 186427= 16 Cyg: hierarchical triple system; 16 Cyg A 
(the companion of the planet host, at a projected separation 
of 850 AU) was shown to be a binary with a projected sep- 
aration of about 70 AU (Patience et al. 2002 1. The masses 
of 16 Cyg A and its companion are about 1.02 and 0.17 Mq 
respectively. 

- HD 189733: this is the only host of a transiting planet in a 
binary system. The plane of the planet orbit is thus known. 
According to Bakos et al. (i2006aj . from the available radial 
velocity and position data of the binary system, there is ev- 
idence that the binary and planetary orbit are not coplanar. 
Therefore, this object is a candidate for the occurrence of 
Kozai migration. Characterization of the binary orbit will al- 
low to test this scenario. 



C.1. Objects not included in B06 catalog 

The criteria for the inclusion of planets in the 'Catalog of 
Nearby Exoplanets' are different to other compilations such as 
the Extrasolar Planet Encyclopedia maintained by J. Schneider^. 
Three stars classified as planet hosts by the Extrasolar Planet 
Encyclopedia and not included in the Catalog of Nearby 
Exoplanets are members of multiple systems: 

- HD 188753 A was announced to host a hot Jupiter by 
Konacki ( 2005b i. This star is a member of a hierarchical 
triple. The close pair HD 188753BC reaches at periastron a 
distance of just 6 AU from HD 188753 A. If confirmed, this 
candidate would be the planet in the closest multiple system 
found up to now. 

- HD 196885: a planet candidate was reported in 2004 at 
http ://exoplanets .org/esp/hd 1 968 85/hd 196885. shtml , with- 
out successive publication of further details. The star was 
found to be a tight binary (0.7 arcsec, 25 AU projected sepa- 
ration) by Chauvin et al. ( I2006> . 

- HD 219449 is a giant star. It was included in the study of R06 
that showed that the system is a hierarchical triple. A close 
pair, HD 219430 (separation 18 AU) is at 2250 AU projected 
separation from the primary (the planet host candidate). 

C.2. Additional unconfirmed binaries 

A few additional planet hosts not included in TableQlhave com- 
panion candidates whose identification or physical association 
needs confirmation: These are: 

- HD 169830: a candidate companion at 1 1 arcsec was shown 
by R06 to have photometric distance compatible with physi- 
cal association with the planet host. However, the low proper 
motion of the primary does not allow confirmation. Further 
observations (radial velocity, spectral classification) are re- 
quired. 

- HD 217107: it is listed in WDS as a binary with a com- 
panion with a projected separation of 0.3-0.5 arcsec=6-10 
AU. The magnitude difference is not provided. Accor ding to 
R06, this detection needs confirmation. Chauvin et al. (I2006I I 
reported null results of their adaptive optics search using 
NACO. Their limit on hKs at 0.3 arcsec is about 7.5 mag, 
corresponding to a mass of about 0.1 M©. Vogt et al. (2005J 
also report non detection of stellar companions at separa- 
tions larger than 0.1 arcsec. The Hipparcos catalog does not 
include evidence for astrometric motion, expected for a stel- 
lar companion at 6 AU. As noted by R06, the presence of a 
stellar companion at 6 AU would not be dynamically com- 
patible with that of the outer planet {a ~ A AU). One pos- 
sible explanation is that the system is seen nearly pole-on 
and that the stellar companion and the outer 'planet' are 
the same object. However, this possibility is not convinc- 
ing considering that the second epoch observation listed in 
WDS is 1997, fairly close to the periastron of the RV orbit 
(1998.7 ± 0.3). Detection should have been easier in the past 
few years, in contrast to the null detections by deep adaptive- 
optics searches. 

- HD 168443: Reffert & Quirrenbach ( l2006l derived an astro- 
metric mass of 36 ± 12 My for the outer 'super-planet' HD 
168443c. As for HD 38529, there is then ambiguity about 
the classification of the object. 

' www.exoplanet.eu 
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- HD 111232 and HD 150706: astrometric acceleration de- 
tected by Hipparcos (see R06). 

- HD 52265, HD 121504, HD 141937, HD 154857, HD 
162020, HD 179949 and HD 183263: first epoch observa- 
tions by Chauvin et al. (2006 1 revealed companion candi- 
dates. Follow-up observations are required to establish their 
physical association with the planet hosts. 
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Table 1. Properties of planets in multiple systems, of their host stars and their companions: planet msin/, period, semimajor axis, 
eccentricity, radial velocity semi-amplitude, metallicity and mass of the planet host (from B06); projected separation, semimajor 
axis and eccentricity (when available) and mass of the companions, and critical semiaxis for dynamical stability of planets. The 
asterisks in the last column mark systems discussed individually in AppendixlO 
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Binary references: HD 142: R06; HD 1237: Chauvin et al. f2006); HD 3651: Mugrauer et al. K006b'); HD 9826: Lowrance et al. (2002 l, E04; 
HD 11964: R06; HD 13445: Lagrange et al. (2006 1, Mugrauer & Neuhauser (2005 I; HD 16141: Mugrauer et al. (2005 1; HD 19994: Hale ( 1994), 
E04; HD 20782: this paper; HD 27442: Chauvin et al. (2006(, R06; HD 38529: R06; HD 40979: E04; HD 41404: E04; HD 46375: Mugrauer et 
al. ( 2006a 1; HD 75289: Mugrauer et al. (2004b 1; HD 75732: E04, Mugrauer et al. ( 2006a 1; HD 80606: E04; HD 89744: Mugrauer et al. (2004^; 
HD 99492: R06; HD 109749: this paper; HD 1 14729: Mugrauer et al. ( 2005 1; HD 1 14762: Patience et al. ( 2002 1, E04; HD 120136: Hale ( 199§, 
E04; HD 142022: Eggenberger et al. (2006b(; HD 147513: Mayor et al. 1 2004.1 . R06, Silvestri et al. ( 2001 1 ; HD 178911: Tokovinin et al. l 2000l . 
E04; HD 186427: Hauser & Marcy ( 1999 1, Patience et al. ( 2002 1, E04; HD 188015: R06; HD 189733: Bakos et al. j^006a 1; HD 190360: E04; HD 
195019: E04; HD 196050: Mugrauer et al. >2005t : HD 213240: Mugrauer et al. <2005l ; HD 222404: Hatzes et al. i2003l : HD 222582: R06. 

" The secondary is a white dwarf; the original mass was larger and the separation smaller. See App.|B]for details. 
Hierarchical triple system, the sum of the masses of the two companions is listed here. See App.|C|for details. 



